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We study the motion of higher dimensional fermions in a non-singular 6D brane background with 
an increasing warp factor. This background acts as a potential well trapping fermions and fields 
of other spins near a 3+1 dimensional brane. By adjusting the shape of this potential well it is 
possible to obtain three normalizable zero mass modes giving a possible higher dimensional solution 
to the fermion generation puzzle. The three different zero mass modes correspond to the different 
angular momentum eigenvalues for rotations around the brane. This bulk angular momentum acts 
as the family or generation number. The three normalizable zero modes have different profiles with 
respect to the bulk, thus by coupling the higher dimensional fermion field to a higher dimensional 
scalar field it is possible to generate both a realistic mass hierarchy and realistic mixings between 
the different families. 

PACS numbers: ll.lO.Kk, 04.50.+h, 11.25.Mj 



I. INTRODUCTION 

In the standard model of particle physics [1] there are open questions which have not yet found an answer. Chief 
among these is the fermion family or generation puzzle as to why the first generation of quarks and Icptons (up quark, 
down quark, electron and electron neutrino) are replicated in two families or generations of increasing mass (the 
second generation consisting of charm quark, strange quark, muon and muon neutrino; the third generation consisting 
of top quark , bottom quark, tau and tau neutrino) In addition to explaining why there are heavier copies of the 
first generation of fermions one would like to explain the mass hierarchy of the generations and the mixings between 
the generations characterized by the CKM (Cabbibo-Kobayashi-Maskawa) matrix. Several ideas have been suggested 
such as a horizontal family symmetry [2] . 

Recently theories with extra dimensions have been used in a novel way to try and explain some of the open questions 
in particle physics and cosmology. In [3] [4] [5] the hierarchy problem (i.e. why the gravitational interaction is many 
orders of magnitude weaker than the strong and electroweak interactions of particle physics) was addressed using 
large or infinite extra dimensions. Early versions of these extra dimensional models were investigated by several 
researchers [()] [7] [>^] ['■)]. In contrast to extra dimensions in the usual Kaluza-Klcin picture, in the models with 
large or infinite extra dimensions gravity acts in all the spacetime dimensions, while the other particles and fields are 
confined, up to some energy scale, to a 3 + 1 dimensional brane. These recent extra dimensional models have also been 
applied to answer other questions of particle physics and cosmology. Non-supersymmetric string models [10] have 
been constructed which have been able to reproduce the standard model particles from intersecting D5-branes. More 
phenomenological brane models [11] [12] [13] [14] have been constructed to explain the hierarchy of masses and/or the 
CKM elements of the fermions. In [15] [16] brane world models were used to explain dark energy and dark matter. 
General studies of higher dimensional cosmologies can be found in [17] [18]. 

In this paper we attempt to give a toy model for the generation problem using a brane world model in 6D. We 
obtain three 4D fermion families from the zero modes of a single 6D spinor field. A mass hierarchy and mixings 
between the three zero modes are obtained by introducing a Yukawa type interaction between the 6D spinor field and 
a 6D scalar field. This gives a common origin (i.e. the higher dimensional Yukawa interaction) for both the mass 
spectrum and the mixings of the fundamental fermions. 

II. 6D GRAVITATIONAL BACKGROUND 

In [13] [19] [20] [21] a 6D brane world model was investigated which gave universal gravitational trapping of fields 
of spins 0, ^, 1, 2, to the brane. The system considered was 6D gravity with a cosmological constant and some matter 



'Electronic address: sma34@csufresno.edu 
tElcctronic address: dougs@csufresno.edu 



2 



field energy-momentum. The action for tiiis system was 



S = J (fxy^- 



-'9 



(1) 



where \J—^g is the determinant, M is the fundamental scale, is the scalar curvature, A is the cosmological constant 
and is the Lagrangian of the matter fields. All of these quantities are six dimensional. The ansatz for the 6D 
metric was taken as 

ds^ = (j>^{r)r]af3ix'')dx°'dx>^ - X{r){dr^ + r^dO^) , (2) 

where the Greek indices a, /3, ... = 0, 1, 2, 3 refer to 4-dimensional coordinates. The metric of ordinary 4-space, r]ap(x'^), 
has the signature (+, — , — , — ). The functions (f>{r) and A(r) depend only on the extra radial coordinate, r, and thus 
are cylindrically symmetric in the transverse polar coordinates (O<r<oo,O<0< 2tt). The ansatz for the 
energy-momentum tensor of the matter fields was taken to have the form 

T'mi/ = -9fj^uF{r), T,j = -gijK{r), 7^^ = . (3) 

Other than satisfying energy-momentum conservation i.e. 

V^Tab = -^dAi^gT^^) + rg^T^^ = K' + 4^iK-F)^0, (4) 



V -'9 

the energy-momentum tensor was unrestricted, although it was desirable for it to satisfy physical requirements such 
as being everywhere finite, and being peaked near the brane. In (4) and in the rest of the paper a prime indicates 
derivative with respect to r 

In [1.3] [20] it was found that the above system had the following non-singular solution 

, , c^' + ar'' 

where a, 6, c are constants and a > 1. All other ansatz functions were given in terms of 4>{r). 

Kr) = ^= p^a l)bc \ , (6) 

where p is an integration constant with units of length, which was related to the constants a and b by 

p^K _ b 
10M4 ^ a- 1 



(7) 



This solution in terms of (f){r) and A(r) represents a non-singular, thick brane. The brane thickness is proportional 
to c. Recently the matching conditions for a general thick brane of codimension 2 was given [22]. This provides a 
general framework in which to study gravitational phenomena and particle trapping in codimension 2 braneworlds. 
The source functions are also determined by 

^("^ = 20^+40' = + 

where /i and /2 are constants given by 

/i = -|^a, /2 = ^(a-fl), (9) 

In [13] and [20] and it was shown that for 6 = 2 this solution gave a universal, gravitational trapping for fields with 
spins 0, |, 1, 2 within the brane width (« c) of r = 0. 

For & > 2 (6) shows that the scale factor for the extra dimensions, A(r) = at r = 0, raising the possibility of 
having a singularity on the brane and making the solutions with b > 2 unphysical. However, by looking at invariants 
such as the Ricci scalar, R, one finds that they are non-singular at r = 0. This indicates that the zero of A(r) at r = 
for b > 2 solutions is not a physical singularity. The Ricci scalar for the above solution is 

26(-5c2'' + Aac^'' + lOcV'' - 22acV' + lOa^cV'' + iar^'' ~ ba^r^'') 

{a-l)p^{c'> + arby ' ^^^^ 
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It is easy to see that this is finite at r = 0. Other invariants such as the fully contracted Riemann tensor, 
Rab '"^Rcd or the square of the Ricci tensor, RabR^^ also turn out to be finite at r = 0. Thus we take 
the zero in the scale factor A(r) to be a coordinate rather than physical singularity. Note that A(r) goes to zero both 
at r = and r = oo so that the metric (2) essentially becomes 4D at these locations. Thus the solution given in (5) 
is like the 2 brane model of [4] where two 4D branes sandwich the higher dimensional, bulk spacetimc. 

The weak point in the above is that the source ansatz functions have no clear physical interprcation. It would be 
desirable to show that an cnergy-mometum tensor of the general form given by F(r),K{r) could arise from some 
realistic source such as a scalar field. Because (j>{r) increases as r ^ oo (8) (9) indicate that F{r),K{r) have their 
maximum near r = and decrease as — ^- oo. This behavior is similar to soliton solutions of classical field theory. Thus 
F{r),K(r) might be considered as modeling some solitonic field configuration which forms the brane. As specific 
examples references [23] [24] investigate solitonic scalar field configurations which form branes. Another possible 
physical interpretation of the source ansatz functions can be given by transforming the metric given by (2) (5) (6) via 
the transformation r = ctan^/''(/3/2) and setting c^/4 = lOAf^/A = p^{a — l)/h. With this the metric takes the form 

ds^ = <l)\(3)r]^f3{x'')dx"dx''^ -jidp"^ + —sm^ I3d9^j , (11) 

where = ^{{a + 1) + (1 — a)cos/3). When 6 = 2 one can see (after renaming the angles in a standard way as 
/3 — > and O^if) that the geometry of the extra dimensions is that of a sphere. However when b > 2 one has zeroes 
in the 2D scale function, A(r), at r = and r = oo and a conical deficit angle of 6 = 2tt — Im. This is essentially 
the "football" -shaped geometry for the extra dimensions considered in [25] (see in particular equation (24) of [25]). 
The 6D solution in [25] involved realistic sources: an electromagnetic field in the form of a magnetic flux and a bulk 
cosmological constant. Thus the source ansatz functions, F(r),K(r) may be acting as effective magnetic flux plus 
bulk cosmological constant. 

III. THREE FAMILIES FROM THREE ZERO MODES 

We now study the motion of a 6D fcrmion field in the gravitational background given by (2). We make the 
identification that different fcrmion zero-mass modes (i.e. solutions for which the 4D part of the fcrmion wavefunction 
satisfies ^^d^il){xy) = 0) correspond to different families. Thus we want to see if it is possible to obtain three zero 
modes. The picture that we present here is a toy model in the sense that the effective 4D fcrmion fields do not 
have the full SU{2>) x SU{2) x U{\) charges of the real Standard Model fields. As will be shown below the fermions 
considered here carry only a U (1) charge associated with the rotational symmetry around the brane i.e. the symmetry 
associated with the extra dimensional variable 9. This is the same as the first example in [12], with the U{\) quantum 
number being associated with family number. Reference [12] gives other more complex and realistic examples where 
the fermions carry U{\) x C/(l) or SU{2) x U{1) charges. In the U{1) x U{\) example one of the t/(l)'s is associated 
with a horizontal, family symmetry [2] which distinguishes between fermions of different generations, while the other 
U{1) is then an ordinary gauge symmetry. In the present case the single UiV) is associated with the family symmetry. 
Other authors [2(i] also have studied 6D brane models with more realistic standard model charges. 

One can gauge [12] any of the above examples using the Kaluza-Klein approach. For example, with the simple U{\) 
model considered here the metric in (2) has a Killing vector, dg, which via the standard Kaluza-Klein mechanism 
implies an associated gauge boson in the effective 4D theory. The gauge field arises from the off-diagonal components 
of the higher dimensional metric as 

ds^ = (j>^{r)i]c,p{x'')dx°'dx^ - X{r){dr'^ + (rd9 + A^dx^^f) . (12) 

This is analogous to the horizontal gauge boson of the family symmetry models [2]. In [27] it was shown that the 
zero mode of the gauge field was localized to the brane. Since horizontal gauge bosons have an experimentally 
fixed lower mass limit a more realistic model would need to have some symmetry breaking, Higgs mechanism in order 
to give Afj^ an acceptably large mass. 

From the previous section we found that in 6D one has the freedom to let the exponent b take values other than 
6 = 2. The motion of fermions in the 6D brane solution of (2) for the 6 = 2 case was studied in [21]. It was found 
that in this case only one zero mode occurred, and thus only one family. Therefore we want to consider the 6 > 2 case 
and show that for a certain range of 6 it is possible to get three zero modes. The constant 6 controls the steepness of 
the scale functions (p{r) and A(r). Therefore it is reasonable that larger 6 should give a stronger confinement of the 
fermions to the vicinity of the brane at r = 0, and thus to a larger number of zero modes. 

The 6D action and resulting equations of motion for a spinor field are 

5* = / d^xy/^^iT'^DA'i , F-^Da* = F^D^^- + rOr'i/ + T^Dg'if = . (13) 
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In the above = ^^7^ Stre the 6D curved spacetime gamma matrices, and /i^ are the sechsbiens defined via 
gAB — f^A^BVAB- III order to evahiate the 6D Dirac equation in (13) we need to calculated the spin connections 



(14) 



The non-zero spin connections are 



With tlrcsc one can explicitly calculate the various covariant derivatives in (13) 

1 



1 



(15) 



(16) 



We now assume that the 6D fermion spinor can be decomposed as \l/(x^) = ip{Xf_i) (8) C{r,9) into 4D and 2D parts. 
We are interested in the zero-mass modes so the 4D fermion part is taken to satisfy 'Jfj,d'^tjj{x^) = 0. The 2D spinor 
can be expanded as 



We take the gamma matrices of the extra space as in [12] [13] 



iW 



(17) 



(18) 



Combining equations (13) - (18) we arrive at the following equations for fi{r) and gi{r) 



b' 1 dr (VW) I 



9^-1-2 



y 1 dr {VW) I 
i> 2 y/rcj)' r 



2 ^/W 

The solutions for /; (?') and gi (r) are 

fi{r) = a,(/)(r)-2(r-0'(r))-3r' , gi{r) = hq^ir)-^ (rcj)' (r))- 



(19) 



(20) 



Because of the I dependence in both fi{r) and gi(r), different I values give fermion fields with different profiles in 
the bulk. This will result in different masses and mixings for different I values when the fermion field is coupled to 
a scalar field in the next section. One criteria for the trapping of the fermion field is that it should be normalizable 
with respect to the extra dimensions r, 9 



1 = / ^^C{r,d)C{r,0)drdd = 



2tt 



d6 dr0VV^(«^"+^'/r-"") = 2V / dr (afr'' + bfr-'') [21) 



From (6) one finds \f^~fr = (a — 1)262027-2 ^(c** -I- r^) ^ . Thus in order for (21) to be normalizable and for the 
particular fermion i-mode to be trapped, we want the integral 



27rpS/(a- l)6c' 



-dr 



(22) 



to be finite. If (22) diverges the particular Z-mode will not be trapped. This requirement that (22) be finite leads to 
restrictions on h for particular values of I. Evaluating the integral (22) gives 



27rV^c=^^' sec ( -j- 



2/7r\ /a-1 



if 6>4|Z|, 



(23) 



and (22) diverges if 6 < 4|/|. Thus in order to have three normalizable i-modes we require that 4 < 6 < 8. Under these 
conditions the / — 0, and \l\ ~ 1 modes are normalized and trapped, while jZj > 2 modes are not. Since the integrand 
in (22) is positive definite and only has possible divergences at r = and r = 00 one can come to this conclusion by 
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investigating the r ^ and r ^ oo behavior of this integrand. For |/| = 1 one finds that for Z = +1 the integrand 
behaves as J?™q ~ , ~ r^~2 ; for / ~ —1, it behaves as 5,™q ~ r-~^+^, J?™^^ ~ r^^^^. The r ^ limit 

of I = +1 and r — > oo limit of Z = —1 give convergent results. On the other hand the r ^ limit of ^ = —1 and 
r ^ oo limit of I = +1 give convergent results only if 6 > 4. One can see the for 5 > 4 the I — mode is normalized. 
For |/| ~ 2 one finds that for / = +2 the integrand behaves as J,*™o — "r^^^ , J.™^^ ~ r^~2 ; for Z = —2, it behaves as 
~ r-5+^ , 5?™^ ~ r-5~5. The I = +2 integral diverges at r ^ oo and the I = -2 diverges at r ^ if 6 < 8. This 
analysis again shows that one has three normalizable modes (i.e. three fermion families) when the A < b < 8. 

One could also consider using the criteria for trapping that the fermion action be finite when integrated over the 
extra dimensions. 




5** = / d^x^/^^iT'^DA'f = 2^ / dr-\ —iafr^^+bfr-'^^) . / d'^x^f^T^i^i-i" d^^j (24) 



The fermions are trapped if the integral over r in the last expression is convergent. This integral is almost the same 
as the last integral in (21). It differs only by a factor of 1/0 which comes from the sechsbien that modifies the 
gamma matrices, 7"^. The explicit expression for 4>~^ \JWh' be read off from (6) and (5). From this one finds that 
(j>^^ y/ <f>' /r (X r^^^(c'' + ar^)^^. The only change with respect to the normalization condition (21) is that r'' ar'' in 
the denominator. Thus the integral of the action over the extra coordinates will have the same convergence properties 
as the normalization condition (21) , thus giving the same conclusion that three zero-mass modes will be trapped if 
4 < 6 < 8. 

In [21] only the 6 = 2 case in (5) was considered and only one zero-mass mode occurred. Thus the existence of 
three zero modes is the result of allowing the exponent in (5) to take values b > 2. In [2S] it was shown that the 
solution of (5) (6) (8) could be generalized to spacetimes of dimension greater than 6D. For these greater than 6D 
spacetimes the exponent, b, in (5) was not free, but fixed to 6 = 2. This would seem to imply that only in 6D can one 
have more than one fermion generation for the background solution given by (5) (6) (8). However in the case where 
spacetime greater than 6D one could consider taking the higher generations as non-zero mass modes. Also one could 
try to generalize the other 6D brane solution given in [21] to spacetime dimensions greater than 6D. 

In discussing the masses and mixing between the different families (i.e. different /) we will need the normalization 
relationship between ai and 6;. From (21) (22) and (23) we find 



This condition allows us to write bi in terms of ai or visa versa. 



IV. MIXINGS AND MASSES 



By adjusting the exponent b in our gravitational background solution we have three zero mass modes which can 
be taken as a toy model for three generations of fermions. There are two problems: first there is no mixing between 
the different generations due to the orthogonality of the angular parts of the higher dimensional wave functions. 
Overlap integrals like J^^ Q(^rndrdO, which characterize the mixing between the different states, vanish since 

/„ e-^'^e^'^d^ = if / 7^ m. Second, aU the states are massless, whereas the fermions of the real world have 
masses that increase with each succeeding family. Following [ ! 2] we address both of these issues by introducing a 
coupling between the 6D fermions and a 6D scalar field of the form Hp{x^)^i{x^)'ifi'{x'^). This adds to the action 
a scalar-fermion interaction of the form 



Ssf=fj d^xdrde^/^Hp^i^ir , (26) 

/ is a constant which gives the scalar-fermion coupling strength. 
We now take the scalar field to be of the form 

Hpix^) = Hp{r)e'P' (27) 

i.e. only depending on the bulk coordinates r,9, but not on the brane coordinates x''. In [12] the same form as in 
(27) was taken for the scalar field, but certain simplifying assumptions were made about the form of H(r) - it was 
assumed to be either a constant or a delta-function. In [14] other forms for the scalar field profile were used. In 
the following we will determine the form of H(r) by studying the field equations for a scalar field in the background 
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provided by (5) (6). Note that in the form (27) the scalar field is only dynamical with respect to the extra dimensions, 
r, 9, but not with respect to the brane spacetime dimensions, x^. Thus one has a scalar field mechanism for fcrmion 
mass generation without a dynamical 4D scalar particle. 
Substituting (27) into (26) we find 

S.f^Uw [ d^xM^ni'i'ixn where Uw ^ f f drd9^/^Hp{r)e:^P~'+''^'0{r)Ci'{r) . (28) 



Uw will be non-zero when p — I + I' = 0. When I = /' this gives a mass term and when I ^ I' this gives a mixing term 
between the I and I' modes. 

To get explicit results for Uw one needs an explicit form for Hp{r). This is done by solving the field equations for 
a test scalar field in the background given by 4){r) and A(r) for the different p— modes. The equation for the scalar 
field in the background given by (5) (6) is 

-^dA (V^g^'^dBHpix^)) = . (29) 



Inserting (27) into (29) we get 



From (5) we find that 

40'(r) _ 4(a- 1)^(0' 



(r) 



(31) 



Inserting this expression back into (30) we were not able to find a closed form solution, but were only able to solve it 
numerically. From the form of Ac/)' (r) / (f){r) in (31) one can sec that when 4 < 6 < 8 that the —p^/r^ and l/r terms 
dominate in the limits r —^ and r —^ oo. If one drops the 40'(r)/(/)(r) term then one finds that the asymptotic 
(r — > and r oo) solutions to (30) are 

Hp{r) = A+pr^P^ or A.pr'^P^ for p^O (32) 
Ho{r)=Ao or Boln(r) for p = 0, (33) 

where A±p, Aq, and Bq constants. These asymptotic solutions gave a fair representation to the numerical solution 
even for intermediate values of r. The p = solutions can be written in combined form as Ho(r) = Bq ln(r/co) where 
Aq = —Bq In(co). This form will be used the next subsection to give masses to the three zero modes. The singularities 
in HQ{r), at r = and r = cx), are not a problem since the combination of the fermion "wave function" and the metric 
ansatz functions go to zero fast enough at r = and r = oo to negate these singularities in the Yukawa coupling 
integral (28). 



A. Masses 



From (28) one can sees that the mass terms are those for which I — I' and thus we want to consider couplings to 
the p = scalar mode in (32). We will use the combined form of the two solutions namely Ho{r) = Bq ln(r/co). With 
this (28) becomes 



r I .oo In f A 

■ j drde^gHo{rMr)Ci{r) = 2/i?oVV(« " I J> + r" 



mi^Uu = f I drde^-^gH^{T)Q{r)Ci{r) = 2fBonp'J{a - l)6c& / ^ , [afr'' + 6f r"^' | 

^ Jo c + r 

/^o(ln(-)+ Vtanf — )) where X, = - 1 . (34) 



Co / b \ b 



In arriving at the final line in (34) where have used (25) to replace bi in terms of a/. Looking at only the tan(2/7r/6) 
term and taking fBoKiir/b > gives a hierarchy of masses of m_i < tuq < m+i. However to_i < and niQ = 
which is phcnomenologically wrong. Taking ln(c/co) as positive (i.e. c > cq) and sufficently large can shift the 
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mass spectrum so that all masses are positive and while still maintaining the hierarchy m_i < mo < Because 
our fermions only carry a U{1) charge the above hierarchy is a toy model. Here, as an example, we take the three 
fermions as the "down" sector of quark generations where / = — 1 is the down quark, Z = is the strange quark, and 
I = +1 is the bottom quark. Taking the masses of the down, strange and bottom quarks as m_i = = 4 MeV , 
Too = ms = 100 MeV and m+i = mb = 4400 MeV we find from (34) 

nid _ m_i _ ^ / -A-_i7rtan(27r/&) ^ \ ^ to^ ^ m+i ^ ^^^^ ^ / K+m tB.n{2n/b) ^ \ ^^^^ 



rris niQ \ 51n(c/co) / ' rUs tjiq \ &ln(c/co) 

Solving these equations for ai and a_i (which arc embedded in the definition of K+i and A'_i) gives 




43a;)cos(yj , a_i = (1 + 0.96a;) cos ) , (36) 



where x = J'^^l^(^2n'/b) ^-"^"-^ ^ \l TJ^- ^'^'^ "^-i chosen as in (36) then the mass ratios in (35) are 

obtained. 



B. Mixings 

A similar analysis can be carried out with the mixings between the different "families" characterized by different 
I number. The mixings are delineated by J/q,!? C^i.Oj C^i.-ij and C/-i,o- In the case of mixings the scalar 

field must have a non-zero angular eigenvalue i.e. Hp{r,9) = Hp{r)e'^P^ with p 0) which satisfies p — I + I' = 0. 
Thus for ?7_i,o and ?7o,i one needs p = —1; for Uifi and J7o,-i one needs p = I; for C/i,-i one needs p = 2; for U-i^i 
one needs p = —2. We will require the following relationship between the mixings: C/o,i = C^i,o, C^i.-i = C^-i.i and 
C/o,-i = U-i,Q. This in turn implies that Hp{r) should depend only on \p\ (e.g. Hi{r) — _ff_i(r)). Looking at the 
first line in (32) this means that we can take cither the first solution or the second but not the sum in general unless 
A^i = A-i. In what follows we will take Hp — A+p?*'^'. The conclusions are not qualitatively different if we make 
the other choice Hp = A_pr~IPl. With this we find 



oo 



f 2 







l + l' 



bibvr 



-i-V 



Uw = / / drde,/^gHp{r)Ci{r)(:v{r) = 2V./\/(a - l)bc^ / Hp{r)-j^^ ajai.r 
where in (37) we have carried out the dO integration, and the condition p — I + I' holds. Inserting Hp{r) = A+pV 



(37) 



(37) and assuming all ai and hi are purely real we get 

rl+P-i r 



2TTp'fd{a-\)hcbA+p 



_|_ 



2TT'p\f 



a-1 



aiai'c''^'' sec 



{l + l' + \p\)7r 



+ bibi'C 



sec 



{l + l'-\p\)^ 



(38) 



Now the four cases (/ = 0, /' = -1), (/ = -1, I' = 0), (/ = 0, /' = 1) and {I = 1, /' = 0) involve \p\ = 1 and from (38) 
yield 



For the two cases / = 1, I' 



Uo,-i = t/-i,o = 
C^o.i 

-1, / = -1, I' = 
Ui,-i - C/-i,i 



fDA+ic^ ( aoa-i 



1 one has |p| = 2 
fDA+2c'' 



bob-i sec 
aoa_i + 



27r 

T 

bob-i 



(aia_i + bib^i) sec 



(39) 
(40) 

(41) 



where as in the previous subsection = ^Jj^^a y ^^-ri • We now show that parameters (i.e. a, b, c, oq) can be chosen so 
the ratios of the above mixings match the ratios of the CKM mixing matrix elements. If this can be done then taking 
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a±i as in (36) will yield the correct mass ratios. From [29] we take (7o,-i = Kts ~ 0.224, C/o,i = Vcb ~ 0.040 and 
U-i^i = Vub ~ 0.0036. Again note that since in our model we only have one flavor in each family (here taken as the 
"down" flavor or sector) these associations between Uij and Vij arc to be taken as representing generic inter-family 
mixing. Now combining (39) (40) (41) gives 



Uq-i _ apa-i + sec (^) bph-i 
Uo,i bobi + sec (^) aoai 



5.6 (42) 



b 

Uo^_, _ A+, + 606-1 sec (^)) 

C/-i,i A+2 (flia-i + 6i&-i)sec (^ 



2^^ 60.5. (43) 



To determine the last ratio one needs to determine the normalization constants. A-^-i and ^+2 of the p = 1 and p ~ 2 
scalar field modes. For this one needs to explicitly evaluate the integral 1 = J \/ —^gH^drdO for p = 1 and p = 2, 
with Hp{r) given by (32). Explicitly the ratio A+1/A+2 is 



A, 



poo .+3 , 

Jo (i+y'>)6 y ay ^^^^ 



The two mass ratios have already been fixed by choosing a±i (this also fixes b±i because of (25)) as in (36). Next 
the mixing ratio, -jj^, hi (42) can be fixed by choosing qq (this also fixes 60 because of (25)). In particular let us 



parameterize as ao = ^rj so that 60 = ^\/2 — ry^ where 77 is arbitrary, and D and c are previously defined. In this 
way one can see that (42) and also (36) are independent of c. Finally one can fix (43) by choosing a in (44) to give 
the ratio A+i/A+2- As an example choose a = 7.5, 6 = 4.041, c = 1.0, D = 1.0 (this can be done by adjusting p), and 
X = 0.001 (this can be done by adjusting cq). In this way one finds ai « 0.1289 and a_i « 0.1263. The associated 
6;'s are 61 0.1234 and 6_i « 0.1262. Using these one find that the mass ratios in (35) are satisfied. Next choosing 
77 w 0.223 so that oq ~ 0.223 and 60 ~ 1.397 one finds that the mixings in (42) and (43) are satisfied. 

Other values of a, 6, c, D and x in this general range worked as well. However in general the various relationships 
worked best when 6 was close to 4. 



V. SUMMARY AND CONCLUSIONS 



In this paper we studied the field equations of fermions in the background of the non-singular, 6D brane solution 
of [13] [19] . By allowing the exponent, 6, in the 4D scale function, (j){r), to take values 6 > 2 wc found that we could 
get multiple zero mass modes which were identified with different fermion generations. In particular for 4 < 6 < 8 we 
obtained three zero mass modes corresponding to different I eigenvalues: I = —1,0, 1. The charge I played the role 
of the family number. When one fixes the value of 6 = 2 as in [19] one has only one zero mode [21]. For 6 > 2 the 
2D scale function, A has a zero both at r = and r = 00. However the scalar invariants such as the Ricci scalar are 
well behaved and non-singular over the entire range of r, indicating these points are coordinate rather than physical 
singularities. 

The masses and mixings between the different generations was given by a common mechanism - the introduction 
of a scalar field with a Yukawa coupling to the fermions. 

An interesting extension of the above scenario is to sec if the scalar field could play a dual role: (i) as the mechanism 
for generating the masses and mixings and (ii) as the matter source for forming the brane. In [23] [24] it was shown 
that a scalar field could be used as a source to construct a thick brane. Thus it might be possible to replace the 
phenomenological matter sources, F{r) and K{r), by a scalar field source. Such a scenario would be more economical 
since the scalar would serve the dual role of forming the brane and giving masses and mixings to the fermions. 
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